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Lateral fluid forces on two-dimensional centrifugal impellers, which whirl on a 
circular orbit in a vaneless diffuser, were reported at the 3rd Workshop in 1984, 
Experiments were further conducted for the cases in which a three-dimensional 
centrifugal impeller, a model of the boiler feed pump, whirls in vaneless and vaned 
diffusers. The influence of the clearance conf iguration between the casing and 
front shroud of the impeller was also investigated. 

The result indicated that the fluid dynamic interaction between the impeller 
and the guide vanes induces quite strong fluctuating fluid forces to the impeller, 
but nevertheless its influence on radial and tangential force components averaged 
over a whirling orbit is relatively small. 

INTRODUCTION 

The lateral fluid forces acting on two-dimensional (2-D) centrifugal impellers, 
which whirl on a circular orbit in a vaneless diffuser, were studied experimentally 
as well as theoretically at the authors’ laboratory, and the results were reported 
at the 1st and 3rd Workshops in 1980 and 1984 [refs. 1 and 2]. The results were 
later summarized and published in references 3 and 4. 

The apparatus for forced whirling motion was then replaced by a newly designed 
one [ref. 5], and the studies were further conducted for a three-dimensional (3-D) 
impeller (model of a prototype boiler feed pump) installed in vaneless and vaned 
diffusers . 

From these experiments the following results are summarized in the present 

paper : 

(1) Comparison of lateral fluid forces on 2-D and 3-D impellers in a vaneless 
diffuser . 

(2) Comparison of lateral fluid forces on the 3-D impeller in vaneless and vaned 

diffusers . 

(3) Comparison of lateral fluid forces on the 3-D impeller for various clearance 

configurations between the casing and frpnt shroud. 

(4) Fluid forces induced by impeller/guidfe vane interaction. 

Similar experimental studies on whirling impeller forces have been conducted at 
California Institute of Technology, Sulzer Brothers, Ltd., and Mitsubishi Heavy 
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Industries, Ltd. A 3-D impeller installed in a volute casing was tested for forced 
orbital motion at the CalTech facility [refs. 6 and 7]. The results indicated that 
the tangential force became excitatory at a low whirl speed ratio, i.e., at high 
supercritical speed, regardless of the pump discharge. 

A model impeller of the boiler feed pump was tested at the Sulzer facility, 
where the test impeller was forced to a unidirectional oscillatory motion [refs. .8 
and 9]. The results were qualitatively similar to those obtained at CalTech, but 
the forces induced by impeller oscillation were considerably larger than those 
obtained at CalTech and at our laboratory. This difference could be attributed to 
the fact that the test impeller at Sulzer was a precise scale model of a real 
configuration, and thus the clearance between the impeller and the surrounding 
casing was kept as small as possible to reduce the disk friction loss. This 
situation resulted in a large restrictive force for any impeller motion in the 
casing. 

Mitsubishi built a new whirling test facility for circular orbital motion and 
tested a 2-D impeller in a volute casing [ref. 10]. The results indicated a 
quantitative tendency similar to those of CalTech and ours, and tangential force 
became excitatory at a low whirl speed ratio regardless of pump discharge. 

The 2-D and 3-D impellers tested in our laboratory indicated less instability. 
The whirl excitation was observed only when the impeller operated at partial 
discharge [refs. 4 and 5]. 

The test data obtained by the above mentioned four facilities are compared 
briefly in Appendix A. A qualitative explanation on the nature of fluid forces on a 
whirling impeller in a vaneless diffuser is given in Appendix B. 

TEST FACILITY 

The test facility is virtually a single-stage, vertical shaft, centrifugal 
pump. The overall test setup, which has been fundamentally unchanged since the 
start of this study, is illustrated in figure 1. The new apparatus for forced 
whirling motion is illustrated in figures 2 and 3. The section views of figures 4 
and 5 describe how the test impeller is installed in the casing for the 2-D and 3-D 
impellers. The experiments on 2-D impellers were done principally to compare them 
with theory, while those on the 3-D impeller were done to obtain impeller data for 
practical applications. The specifications of the test impellers and diffusers are 
listed in tables 1 and 2. 

In the previous test facility, the main pump shaft was driven by a constant 
speed induction motor (530 rpm), while the whirl shaft was driven by a variable 
speed DC motor, so that the ratio of whirl speed Q to rotational speed co, Q/w, 
could be arbitrarily adjusted and kept constant. However, this made it impossible 
to keep the phase relation between whirl and rotation angles, 9^ and (see 
fig. 6), in a prescribed way. 

To improve this situation, the DC motor for whirling motion was replaced by an 
AC servo motor. After counting signals from encoders attached to the pump shaft 
(1500 ppr ) and whirl shaft (3000 ppr), proper pulse series are fed to the controller 
of the servo motor, which keeps the whirl angle ©]_ in a prescribed relation to the 
rotation angle 02- At Q/co = 1/2, for example, the impeller rotates twice and 
returns exactly to the same position relative to the guide vanes while the whirl 
shaft rotates once. By introducing this equipment, the impel ler /guide vane 
interaction, where sensitivity varies according to the relative position between the 
impeller and the guide vanes, can be measured with a better resolution. 
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DATA EVALUATION 


The data acquisition and evaluation were made fundamentally in the same way as 
described in the previous reports. The fluid force F(0 1 ,©2) is divided into 
radial and tangential components, F r and F q, as illustrated in figure 6. The 
average forces on the orbit, F r and Fq, are then normalized by the equation, 

f r = F r /Meco 2 , £q = FQ/Meco 2 (1) 

where reference mass M = ptr r 2 2 b, e = eccentricity, p = fluid density, C 2 = outer 
radius and b = impeller exit width. 

In comparing normalized fluid forces on 2-D and 3-D impellers, the following 
difference must be kept in mind. 

(1) For 2-D impellers, fluid forces acting on impeller blades, excluding all other 
fluid forces on front and back shrouds, impeller hub, etc., are normalized by Mew 2 , 
where b of the reference mass M is the exit width of the impeller blades only. 
Reference mass is therefore the mass of fluid in the volume enveloped by the 
rotating impeller blades. 

(2) For 3-D impellers, fluid forces acting on the whole impeller (blades as well as 
shrouds and hub) are normalized by the quantity Mew 2 , where b is the exit total 
width of the impeller, that is, blade width plus the width of the two shrouds. 

The former definition is convenient for comparing the experimental result with 
the theoretical one. The latter is for the convenience of practical application, 
since the fluid forces on impeller blades and shrouds have the same significance on 
the rotordynamic behavior of the rotating system. 


FLUID FORCES ON 2-D AND 3-D IMPELLERS IN VANELESS DIFFUSERS 

Comparison of the two impellers is dealt with in this chapter. The first is a 

2- D impeller whirling in a vaneless diffuser or in free space as shown in figure 4. 
The other is a 3-D impeller whirling in the casing as shown in figure 5. The 
discharge from the impeller is diffused in the succeeding parallel walled vaneless 
diffuser (11 vanes are removed). The clearance between the casing and front shroud 
of the impeller is sufficiently large (without the spacer), and the effect of the 
fluid dynamic restriction due to the surrounding casing is seemingly weak. 

Figures 7 (a) and (b) show how normalized tangential and radial force com- 
ponents of 2-D and 3-D impellers vary with whirl speed ratio, Q/w, at design 
discharge $/$D = 1*0 and at partial discharge <t>/4>D = 0.2. As seen from these 
figures, lateral fluid forces on the 3-D impeller are slightly larger than those on 
the 2-D impeller. Most of this difference could be attributed to the fact that the 

3- D impeller has an inlet passage that is wider than the exit and that the real 
volume enveloped by the rotating impeller is considerably larger than that of the 
reference mass, irr 2 2 b. In other words, the forces on the 3-D impeller are 
normalized by a relatively smaller quantity compared with the case of the 2-D 
impeller . 

It can be concluded from the results that fluid forces on whirling impellers in 
a vaneless diffuser, reqardless of 2-D or 3-D, are caused mostly by unsteady fluid 
forces on the whirling impeller blades and coincide with those predicted 
quantitatively by theory [refs. 3 and 4]. This is only the case, however, when the 
impeller is installed in the casing with ample clearance. 
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FLUID FORCES ON 3-D IMPELLER IN VANELESS AND VANED DIFFUSERS 


In this chapter the influence of the interaction between the impeller blades 
and guide vanes on the whirling fluid forces is dealt with. A 3-D test impeller was 
installed in the casing as shown in figure 5. The eleven guide vanes attached 
between the two parallel walls (see fig. 10) are removable, and whirling tests were 
conducted for cases with and without guide vanes. The mean radial clearance between 
the impeller outer periphery and the guide vane inner periphery is 5 mm. When the 
impeller whirls with eccentricity e = 1.5 mm, the clearance varies between 3.5 and 
6.5 mm. Every time a trailing edge of an impeller blade passes in the vicinity of a 
leading edge of a guide vane, a fluid dynamic interaction takes place, and the blade 
receives an impulsive force. The impel ler /guide vane interaction can be understood 
as the sum of these individual impulses. The clearance between the casing and front 
shroud is kept wide (without a spacer) throughout this test series. 

Figures 8 (a) and (b) show the comparison at design discharge 4>/4>d = 1.0 and 
at partial discharge = 0.2. It is obvious that the fluid forces on an 
impeller in a vaned diffuser are about 50% larger than those on an impeller in a 
vaneless diffuser. It should be noted that the pressure rise in a vaned diffuser is 
also larger by about 15 to 20% at the same discharge. The general qualitative 
tendencies of both cases are essentially the same. As described in the latter 
chapter, impeller/guide vane interaction generates quite large fluctuating fluid 
forces on the whirling impeller; however, its influence on fluid forces averaged 
over an orbit remains relatively small. 

As seen in figure 8 (b), both tangential and radial components indicate a local 
increase at around Q/co = 0.1 to 0.2 at partial discharge. These phenomena are 
hardly seen near design discharge and are considered as characteristics of partial 
discharge. The real reason is still unknown. The coupling of the whirling motion 
with the rotating stall in the diffuser, which is usually 10 to 20% of rotational 
speed, can be considered as a possible cause. 


INFLUENCE OF CLEARANCE BETWEEN CASING AND FRONT SHROUD 

The influence of clearance between the casing and front shroud of the impeller 
on whirling fluid forces is dealt with in this chapter. Tests were conducted for a 
3-D impeller in three cases: (a) without a spacer (original layout); (b) with a 
spacer. Cl = 2.7 mm, C2 = 3.0 mm (see fig. 5); (c) 24 radial grooves with 10 mm 
width and 5 mm depth are added to the spacer surface facing the front shroud. Other 
conditions are the same as those of case (b). 

The result obtained at design discharge is summarized in figure 9. As the 
clearance between the casing and impeller shroud gets smaller and thus the 
restrictive effect of the fluid contained in the narrow space becomes greater, the 
fluid forces on the front shroud of the whirling impeller increase. This results in 
the increase of the radial component and the decrease of the tangential one. 
Narrower clearance between the impeller and casing leads obviously to the expansion 
of the whirl excitation range at positive whirl. 

The grooves on the spacer surface, which are expected to brake the rotation of 
fluid in the clearance, resulted in the increase of the tangential component in the 
positive direction, thus reducing the range of whirl excitation at positive whirl. 
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IMPELLER/GUIDE VANE INTERACTION 


The fluid forces on the whirling impeller are a function of both the whirl 
angle 0^ and rotation angle 02* As understood from figure 10, which illustrates 
the relative disposition of 6 impeller blades and 11 guide vanes, the fluid force at 
a location of whirl orbit is different according to the relative angular position of 
the impeller to the guide vanes. Since six impeller blades are equally located on 
the periphery, the rotational condition becomes identical after every rotation of 
60°. 

Figure 11 illustrates how fluid force on the whirling impeller varies on the 
orbit because of the impeller/guide vane interaction. The measurement was done at 
Q/o) = 1 and at four discharge conditions: that is, 4>/<j>p = 1.4, 1.0, 0.6, and 0 
(shut off). At Q/co = 1 the trailing edge of a specific impeller blade with the 
narrowest gap passes successively near the leading edges of the 11 guide vanes 
during one rotation. This is the reason why 11 distinct interactions can be 
observed during an orbit. As seen from the figures the interaction becomes 
obviously weaker as the flow rate decreases. At partial discharge the separation in 
the impeller and diffuser causes large random fluctuations of flow, and it applies a 
strong masking effect on the blade interaction. 

Figure 12 illustrates fluid force vectors on the orbit for various whirl speed 
ratios, Q/w = 2/3, 1/2, -2/3, and -1. The figures are for the over-discharge 
condition, ^/^D = l* 4 ' where the interaction takes place most distinctly. The 
frequency of periodicity on the orbit varies with the whirl speed ratio and can be 
explained by tracing the mode of the impeller/guide vane encounter near the minimum 
gap region. 

The impeller/guide vane interaction takes place at every encounter of moving 
blade and stationary vane, however, the resultant fluid force on the whole impeller 
is predominantly induced by the blade interaction to a specific guide vane, which is 
located closest to the least gap point. As illustrated in figure 10, we chose the 
guide vane with the least gap and the blade closest to the vane. The angle a is 
then defined as the angle difference between the trailing edge of the blade and the 
leading edge of the vane. All interaction data at 4>/4>p = are then replotted 
against this angle a for three whirl speed ratios, Q/oo = 1.0, 1/2, and 1/6, as 
shown in figure 13, where the tangential and radial components indicate the 
deviation from the average values on the orbit. This figure suggests that the 
interactions at various whirl speed ratios are almost identical in nature and are 
primarily controlled by the local geometric condition near the minimum gap point. 

It should be noted again that the interaction, though it appears quite 
distinctly, is only a local event and has a relatively small effect on the fluid 
force averaged over the orbit. 


CONCLUSIONS 

Fluid forces on whirling centrifugal impellers are studied experimentally for 
various conf igurations of the impeller and diffuser. The results obtained are 
summarized as follows: 

(1) Comparison of fluid forces on 2-D and 3-D impellers whirling in a vaneless 
diffuser indicated no significant difference. 

(2) Fluid forces on a 3-D impeller whirling in a vaned diffuser were about 50% 
larger than those in a vaneless diffuser, partly due to the interaction between the 
impeller blades and guide vanes and partly to the increased pressure rise. 

(3) The fluid forces acting on the front shroud of the whirling impeller are 
affected by the clearance configuration between the casing and shroud. The smaller 
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the clearance, the bigger the fluid forces on the whirling shroud. Grooves cut on 
the spacer surface facing the shroud do augment the damping effect at positive whirl. 
(4) Impeller/guide vane interaction generates strong fluctuating fluid forces on the 
impeller. The effect of this strong interaction on fluid force characteristics is 
relatively mild as stated above, because high-frequency fluctuations are smoothed 
out during the averaging process. 


F : lateral fluid force 

f : normalized fluid force 

M : reference mass 

: whirl angle 

©2 : rotation angle 


SYMBOLS 

4> : flow coefficient 

Q, co : whirl and rotational frequency 

subscript 

D : design (shock-free entry) 

r, 9 : radial and tangential 


APPENDIX A: COMPARISON OF TEST RESULTS 

The facilities and test results of four research groups, that is, the 
University of Tokyo, California Institute of Technology, Sulzer Brothers, Ltd., and 
Mitsubishi Heavy Industries, Ltd., are introduced briefly as follows: 

Univ. Tokyo The test facility is described in the main text and is shown in figures 
1 through 5. Normalized tangential and radial components, fQ and f r , on the blades 
of a whirling 2-D impeller in a vaneless diffuser are shown in figure A-l, in which 
calculated results are also shown for comparison [ref. 4]. The results of a 3-D 
impeller whirling in a vaned diffuser are shown in figure A-2 for design flow rate 
and shut-off. The thin broken curves in the figure indicate the force 
characteristics obtained by the best fit of mass, damping, and stiffness matrices. 
In comparing the results of the four groups, it should be noted that the 
normalization of tangential and radial forces is done by the same definition; 
however, the sign of the tangential force defined by the Univ. Tokyo group is 
opposite to that of the remaining three groups. 

CalTech The test facility is illustrated in figure A-3, in which a horizontal -shaft 
3-D impeller installed in a volute casing is forced to a circular orbital motion and 
the fluid force on the impeller is measured by a rotating dynamometer. Typical 
results at design flow rate are reproduced in figure A-4 from reference 6. 

Sulzer Figure A-5 illustrates the test facility, in which a horizontal-shaft boiler 
feed pump impeller in a vaned diffuser is forced to a unidirectional motion. 
Impeller force is measured by a rotating dynamometer. The tangential force is 
evaluated indirectly from the reaction force to unidirectional excitation and is 
given in figure A-6 [ref. 8]. 

Mitsubishi Figure A-7 illustrates the test facility, in which a vertical-shaft 2-D 
impeller in a volute casing is forced to a circular orbital motion. The results are 
shown in figure A-8 [ref. 10]. 

The comparison of the above results leads to 

(1) Very roughly speaking, tangential force is nearly linear with respect to whirl 
speed ratio, while radial force is parabolic with its minimum at positive whirl 
speed ratio in all test results. 

(2) The absolute values of fluid forces are different. The force level increases in 
order of Univ. Tokyo, Mitsubishi, CalTech, and Sulzer. The forces on the Sulzer 
impeller are almost 8 times as large as those on the blades of the 2-D impeller 
measured by the Univ. Tokyo. 
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(3) The large difference of measured fluid forces comes mostly from the difference 
of clearance space between the impeller shroud and the surrounding casing. The 
narrower the clearance, the stronger the restrictive effect of the surrounding 
casing to the movement of the impeller, thus resulting in a larger reaction force on 
the impeller. 

(4) The results of CalTech, Sulzer, and Mitsubishi indicate that tangential force 
becomes destabilizing at a smaller positive whirl speed ratio, Q/co < 0.4, regardless 
of pump discharge. The results of the Univ. Tokyo indicate, however, that the 
destabilizing effect occurs only at partial discharge. 

(5) The characteristics of whirling fluid forces can be represented by a linear 
formulation using mass, damping, and stiffness matrices to a satisfactory extent. 
However, peculiar behaviors such as the local increase of fluid forces observed at 
low whirl speed ratio and partial discharge (cf. figure A-2) are all smoothed out 
and neglected, when the linear formulation is adopted. 


APPENDIX B: QUALITATIVE EXPLANATION OF WHIRLING IMPELLER FORCE 

The questions of why whirling fluid forces occur and how they are affected by 
whirl speed are explained here qualitatively by using a very simple model. Let us 
consider the fluid dynamic forces on blades A and B, which are indicated by solid 
straight lines in figure B-l (a). It is assumed that the fluid dynamic character- 
istics of the blades are approximated by those of an isolated airfoil neglecting 
interaction between blades. The figure shows the impeller geometry at the moment 
when absolute center O’, impeller center 0, and two blades are aligned on a 
straight line: that is, on the x-axis in this case. This condition is not a 
singular one. For an impeller with a large number of blades, there are always a 
couple of blades that are at the same disposition as that of the figure. The two 
blades considered in the model symbolize the ones which happen to be in such a 
situation. 

When there is no whirl, the lift on a blade, Lq, is proportional to the square 
of peripheral speed u = Rw, where R is the radius of the blade center, and its 
direction is normal to the relative inlet flow W as seen in the figure. If whirl 
occurs, the impeller center 0 circles along the orbit with velocity e£2. At the 
location shown in the figure the whole impeller moves upward with the same velocity, 
so that the absolute peripheral speed of blade A increases from u to u+cQ, while 
that of B decreases from u to u-eft. Consequently the lift of blades A and B vary 
by ±AL, where AL is estimated by a small perturbation principle as AL = 2L 0 eQ/u. 

The imbalance of lifts on blades A and B results in an impeller force Fl shown 
in figure B-l (b), which is proportional to whirl velocity eQ. A body in an 
accelerative motion receives an inertial reaction from the surrounding fluid 
equivalent to apparent mass times acceleration. Let M a be the apparent mass of 
impeller blades, the inertial force of the whirl motion is M a eQ 2 in the direction 
normal to the orbit. This inertial force is indicated by F M in the figure. The 
total fluid force F caused by the whirl motion is therefore the resultant of two 
force vectors, F M and F L . The radial and tangential components of force F are 
then written by 

F r = M a eQ 2 - k r eQ, F q = k@eQ (Bl) 

where k r and k@ are proportional constants. 
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The force characteristics predicted by the above modelling are illustrated in 
figure B-2 by solid curves. The measured characteristics indicated by thin broken 
curves are generally similar in shape but do not necessarily cross the origin as 
solid curves do. Radial and tangential forces at zero whirl speed ratio correspond 
to direct and cross stiffness, respectively. The above analysis, though it is guite 
simple and qualitative, suggests the mechanism of whirling fluid force and also its 
relation to the whirl speed. 

As seen from figure B-l (a), the absolute center O', which represents the 
center of suction pipe and thus the location of source Q, changes its relative 
position to impeller blades during whirl. This induces a source-oriented 
perturbation on the blades and results in a corresponding unsteady fluid force. 
This effect has only a secondary significance as compared to the primary effect due 
to whirl speed-oriented perturbation, however, it constitutes a part of the 
discrepancy between the model and reality. 
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Table 1 Specifications of 2-D impeller 


Impeller : 

Two-dimensional, closed type 

outer blade diameter; 

350 mm 

inner blade diameter; 

175 mm 

exit blade width; 

35 mm 

exit total width; 

45 mm 

number of blades; 

6 

exit vane angle to radius; 68°(logarith. spiral) 

mass; 

4.1 kg 

Diffuser: 

Parallel walled, vaneless 

exit diameter; 

700 mm 

width; 

70 mm 


Table 2 Specifications of 3-D impeller 


Impeller: Three-dimensional, closed type [ 

outer blade diameter; 

330 mm 

inner blade diameter; 

180 mm 

exit blade width; 

28 mm 

exit total width; 

43.5 mm 

number of blades; 

6° 

exit vane angle to radius; 

65.5 

mass ; 

5.6 kg 

Diffuser: Parallel walled, 

vaned and vaneless 

inlet diameter; 

340 mm 

exit diameter; 

480 mm 

width; 

30 mm 

number of vanes; 

11 



Figure 1. Layout of test facility 





Figure 3. Layout of whirling apparatus 
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Figure 4. Casing configuration of 2-D impeller 



Figure 5. Casing configuration of 3-D impeller 
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Figure 6. Coordinate system and force components 
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Figure 10 • Interaction between 
impeller blades and guide vanes 
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Figure 9. Influence of clearance 
configuration on fluid forces of 
3-D impeller 
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Figure 12. Impeller/guide vane interaction on the orbit 
at <t>/4> D = 1.4 
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Figure 13. Force fluctuation plotted against a: 
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A-3. Test facility of CalTech 
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